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ABSTRACT 

We investigate the extent to which cosmic size magnification may be used to com- 
plement cosmic shear in weak gravitational lcnsing surveys, with a view to obtaining 
high-precision estimates of cosmological parameters. Using simulated galaxy images, 
we find that size estimation can be an excellent complement, finding that unbiased 
estimation of the convergence field is possible with galaxies with angular sizes larger 
than the point-spread function (PSF) and signal-to-noise ratio in excess of 10. The 
statistical power is similar to, but not quite as good as, cosmic shear, and it is subject 
to different systematic effects. Application to ground-based data will be challeng- 
ing, with relatively large empirical corrections required to account for with biases for 
galaxies which are smaller than the PSF, but for space- based data with 0.1 arcsecond 
resolution, the size distribution of galaxies brighter than i ~ 24 is ideal for accurate 
estimation of cosmic size magnification. 
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1 INTRODUCTION 

General Relativity predicts that the path of light is distorted 
when interacting with a gravitational field. This modifica- 
tion of the light paths is called gravitational lensing and 
is a powerful tool for probing the distribution of mass in 
the Universe. The variation of the light path depends on 
the position in the sky of the emitting object, the distance 
from the emitting object to the lens and the observer and 
on the potential along the light path. As the Universe is 
in permanent evolution photons that come from an ear- 
lier stage of the Universe will be deflected in a different 



ing see 


Schneider et al.lll992 


;lNaravan & Bartelmannlll996l; 


iMellierl 


19991; iMunshi et al.l 


2008). Hence by studying the 



lensing observable distortion, a reconstruction of the three- 
dimensional unbiased distribution of matter (both light and 
dark) can be performed. The result is a picture of the Uni- 
verse over time, therefore both the power spectrum and the 
growth of density perturbations with redshift can be in- 
ferred. Dark energy and modifications to Einstein gravity 
act to modify the lensing effect by changing the distance- 
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redshift relation and the growth of density perturbations 
l|Hutererll200"2l ; IMunshi fc Wang||2005n . Therefore lensing ef- 
fects are a source of valuable information for three of the 
important open issues in modern cosmology, namely the dis- 
tribution of dark matter, the properties of dark energy and 
the nature of gravity. 

The amount of gravitational lensing induced on the im- 
age of a background source is characterized by an image 
distortion matrix, which may be written in terms of the 
convergence k, that is the field responsible for the changes 
in the image size and brightness, and the complex shear 7 
that, in the limit of small k, parameterises the distortion 
of the shape. Observationally there are three main effects 
on the background sources: changes in the ellipticity, am- 
plification in terms of flux, and magnification of the size. 
For large densities of (dark) matter the effects are all very 
strong, and multiple images of the background galaxies can 
be produced with large distortion. The studies of this dis- 
tortion have led to local reconstruction of the distributio n 
of non-baryonic matter for many years (|Tvson et al.lll984h . 
When the concentration of matter is not as high, the effects 
are not as obvious, and more difficult to detect. The lensed 
background source is observed with changes in the elliptic- 
ities, flux and size, but the effects are so small that, for a 
single galaxy, it is impossible to know if they are intrinsic 
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properties of the galaxy or they are produced by lensing. 
However in this weak lensing limit on large scales, the dis- 
tortions may be observed statistically using a large sample 
of sources. 

In weak lensing, the most studied effect is the modi- 
fication to the galaxy shape (complex ellipticity), a mea- 
sure of the shear. The shape distortion has the main ad- 
vantage that the intrinsic distribution of galaxy elliptici- 
ties is expected to be random, according to the cosmolog- 
ical principle. Therefore under this assumption the aver- 
age complex ellipticity is zero. Weak lensing effects using 
galaxy ellipticities have been detected by several groups us- 



IWittman et al.l (1200(f); 


Semboloni et al. (l2006l);|Jarvis et al.l 


(2006);lBeniamin et al. 


( 20071) ;ISchrabback et alj(2010h. An 



important effort has been made to include and test many 
possible systematic effects on shape measurement (includ- 
ing the point spread function [PSF], instrumental noise, 
pixelization for example), and there are several algorithms 
that can meas ure shapes w i th va rying degr ees of accuracy 
such as KSB i|Kaiser et ail I1995T), KSB+ ([Hoekstra et all 



ll99ST ) and its variants (TRhodes et all |2000l; iKaiserl \200(j) 
and shapelets ([Bernstein fc Jarvisll2002i ; iRefregier fc BacorJ 
2003; Kuiikcn 2006) amongst others. To test, in a blind way, 
the ability of methods to measure the shapes of galaxies a 
series of simulations have been created: STEP, GREAT08 
and GREAT10, where several methods have been tested 
and compared systematically, an exp lanation of the method s 
and their per fo rman ce are shown inlHevmans et alj ( 2006); 
iMassev et ail l|2007h ; iBridle et al] <|2010l ); iKitching et al l 



( 20121 ) respectively. A no vel Bayesian model fitting approach 



lens&t was presented in iMiller et al.l |2007) , that includes 
PS F and pixelizat i on eff ects, and was successfully applied 
by IKitching et~ail (|2010t ) to the STEP challenge images. 
Moving beyond measuring the shapes of galaxies, the am- 
plification effect is usually quantified usi ng the number of 
coun t s of a given flux (see for ex ample iBroadhurst et al.l 
1 19951 ; iHildebrandt et all 2009 , 2011); the number of observed 
galaxies at a given flux increases due to the weak lensing of 
the foreground galaxies, we do not study this further in this 
paper. 

In contrast to galaxy ellipticity measurement, the size- 
magnification effect has not been studied in detail, possibly 
because the complicating effects of the PSF and pixellisation 
were thought to be too challenging. However, there are two 
reasons for revisiting size magnification as a potential tool 
for cosmology: one is that accurate shear estimation is itself 
very challenging, and size could add useful complementary 
information; the second is that methods devised for elliptic- 
ity estimation must deal with the PSF and pixellisation, and 
as a byproduct provide a size estimate, or a full posterior 
probability distribution for estimated size, which is currently 
ignored or marginalised over. In terms of signal-to-noise of 
shear or convergence estimation, the relative strengths of the 
methods depend on the prior distributions of ellipticity and 
size. The former has an r.m.s. of around 0.3-0.4; for bright 
galaxies, the Sloan Digital Sky Survey (SDSS) found that 
the size distribution of massive galaxies is approximately 
log- normal with ulni? ~ 0.3, and for less massive galax- 
ies cr(ln R) ~ 0.5, where R is the Petrosian half-light radius 
l|Shen et al l 120031). For deeper spac e data the dispersion is 
modestly larger jsimard etaL 2002). Thus one might expect 



a slightly a smaller S/N for lensing measurements based on 
size r ather than ellipticity , but not markedly so. Further- 
more, iBertin fc Lombardil (|2006T ) proposed a method based 
on the tight relation between sizes and the central veloc- 
ity dispersion (the fundame ntal plane relat i on) t o reduce 
the observable size variance. iHuff fc Gravel l|201ll ) applied 
a similar method to 55,000 galaxies of the SSDS catalogue, 
and find consistency with shear using the same sample. Also 
a detection with COSMOS HST surv ey using a revised ver - 
sion of the KSB method is claimed in lSchmidt et al] l|2012i ). 

Revisiting the size magnification in detail is the aim of 
this paper. As we will show, to use size-magnification we 
require i) a large area survey to result in sufficiently high 
number densities of galaxies to overcome the intrinsic scat- 
ter, and ii) a consistently small PSF that does not destroy 
the size information of the observed galaxies. Both of these 
requirements can be met with a wide-area space-based sur- 
vey, although some science may be poss ible from the ground . 
Space-based surveys such as EuclicQ (|Laureiis et al.l [201lh 
should meet these requirements (large samples will be avail- 
able, and the PSF size is smaller than typical galaxies), so 
the size information should be considered as a complemen- 
tary cosmological probe to weak lensing ellipticity measure- 
ments. 

One advantage of using the size information is that the 
magnification is directly related to the convergence field k, 
which is in turn directly related to the weighted surface 
mass density, while 7 is related to the differential surface 
mass density. Those different ra dial dependences can b e very 
useful to lift any degeneracies (R ozo fc SchmidtllioTol ). and 
combining the size magnification information with the shear 
will reduce the uncerta inties on the reconstruction of the 
distribution of matter (|jaml 120021 ; I Rozo fc Schmidtl lioiol ; 
IVallinotto et a.1.1 l201lh . All of the shear estimation meth- 
ods referred to above already estimate the size of galaxies 
when calculating the ellipticities, so we expect to measure 
this additional information for free, given an accurate el- 
lipticity measurement, but the accuracy of size information 
should not be taken for granted: it is important to know the 
uncertainties in size measurement, and how they propagate 
to a convergence field estimation. It is this question of how 
accurately one can measure the sizes of galaxies, that this 
paper will address. 

The paper is organized as follows. First, in section 2 wc 
will present the estimator that we will use, then a brief com- 
ment on the method and the characteristics of the simulated 
images. In section 3 the analysis and results are explained 
and finally we will summarize the conclusions in section 4. 



2 METHOD 

A good algorithm for ellipticity weak lensing analysis must 
be able to measure shapes or sizes with an accuracy of a 
few percent, taking into account any possible systematic er- 
rors as the distortion introduced by the PSF, shot noise or 
pixelization. Another requirement is that it should be fast 
because the statistical analysis will be made on large sam- 
ples. This means that the algorithm development is chal- 
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longing because of the dissonant requirements of both in- 
creased accuracy and increased speed as the required sys- 
tematic level decreases. Several methods have been pro- 
posed and applied to weak lensing surveys, in particular 
several of these methods are decribed in t he challenge re- 
ports of STEP, GREAT08 and GREAT10 dHevmans et all 
l200d : iBridle et al.ll20ld ; iKitching et al.ll2012T i. These blind 
challenges have been critical in demonstrating that methods 
can achieve the required accuracy for upcoming surveys by 
creating simulations with controlled inputs against which 
results can be tested. Here we propose a very similar ap- 
proach as the one presented in the GREAT10 challenge, we 
have used simulated galaxy images with different properties 
to measure the response of the size/convergence measure- 
ment under different conditions (corresponding to changes 
in the PSF, S/N and bulge fraction). 



2.1 Estimator 

If the intrinsic size of a galaxy s 3 is magnified, s = /j,s s , where 
to first order the magnification is related to the convergence 
like n = 1 + k, then « = -4- — 1. We can construct an 
estimator for k in the weak field limit by assuming that the 
mean size value is not modified by lensing, i.e., {s 3 } = (s): 



k = — - I. 



(1) 



This assumption should be valid because («) = 0. We cannot 
obtain accurate estimates of k for a single galaxy, we can see 
that from the definition of the estimator, smaller galaxies 
than the mean will always give a negative k, while larger 
galaxies will produce positive k. What is important is to test 
if our estimator is unbiased over a population to a sufficient 
degree to be useful for real data. 

2.2 lensfit 



Throughout we use lenst it (jMiller et al.ll2007l ; IKitching et al.1 
|200S| ; IMiller et al.ll2012T lto estimate the galaxy size; we use 
this because: 1) it has been shown that lens&t performs well 
on ellipticity measurement; 2) it is a model fitting code 
such that it measures the sizes of galaxies (that are part 
of the model); 3) it allows for the consistent investigation 
of the intrinsic distribution of galaxy sizes through the in- 
clusion of a prior on size, and 4) it includes t he effects of 
PSF a nd pixellisation. lens&t was proposed in IMiller et al.l 
( 20071) and has been proved to be a successful tool for g alaxy 
ellipticity shape measurements (|Kitching et al.l 120081 ). Al- 
though model-fitting is the optimal approach for this type 
of problem if the model used is an accurate representation 
of the data, the main disadvantage is that is usually com- 
putationally demanding to explore a large parameter space. 
lens&t solves this problem by analytically marginalizing over 
some parameters that are not of interest for weak lensing 
ellipticity measurements, such as position, surface bright- 
ness and bulge-fraction. The size reported by lensfit is also 
marginalised over the galaxy ellipticity. 

2.2.1 Sensitivity correction 

IMiller et al.l (|2007T l introduced the shear sensitivity, a factor 
that corrects for the fact that the code measures elliptici- 



ties but that shear (a statistical change in ellipticity) is the 
quantity of interest. A similar correction is required for size 
measurement, whereby we measure the size but it is the con- 
vergence that is the quantity of interest; this correction is 
needed because for a single galaxy the prior information for 
the convergence is not known, and we assume it is zero. With 
a Bayesian method we can estimate the magnitude of this ef- 
fect for each galaxy, a further reason to use a Bayesian model 
fitting code in these investigations. Consider the Bayesian 
estimate of the size of galaxy i and write its dependence on 
k as a Taylor expansion: 



Si — Si -\- K 



dsi 



(2) 



In the simplistic case where the likelihood C(s) is described 
by a Gaussian distribution with variance b 2 , with an ex- 
pected value s, and a prior V(s) that also follows a Gaussian 
distribution centred on s with variance a 2 , then the posterior 
probability will follow a Gaussian distribution with expected 
value: 



(s) 



sb + s d a 
a 2 + b 2 



and variance 



2 7 2 

2 _ a b 
° ~ a 2 + b 2 

Note that here s stands for the fitted model parameter for 
the size explained in Sec. 12. 31 and Sd is the data. These equa- 
tions illustrate that the posterior is driven towards the prior 
in the low S/N limit (b — > 00), and thus requires correction. 
Differentiating the expression Sd = s 3 (l + k) we find that 
the k sensitivity correction is: 



ds 
<Ik 



dsn 



a 2 + b 2 dn a 2 + b 2 
substituting into eq. [2] 



(3) 



a 2 + b 2 



we find the estimator for k will be the same as in eq. Q 
corrected by the sensitivity factor: 



s 



a 2 + b 2 



(4) 



In this work we have used this approximation for simplicity 
but in general a normal distribution should not be assumed. 
A more general estimation of the k correction can be done 
in the same way as with the shear and can be evaluated 
numerically, without the need of using external simulations. 

To calculate the sensitivity correction in the general 
case we consider the response of the posterior to a small 
K, £(s — s s ) 1 — ^ C(s — s s — ks 3 ) and expand it as a Taylor 
series: 



C(s — s s — ns s ) ~ C(s — s s ) — 
We then substitute into 

JsV(s)C(s)ds 



dC 

ds 



JP(s)C(s)ds 



and differentiate to obtain the analytic expression for the 
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k sensitivity (for more details of this applied to elliptic ity 
measurement see lMiller et al.ll2007i ; IKitching et alj|2008h 



dn JV{s)£(s)ds 



(5) 



If the prior and likelihood are described by a normal dis- 
tribution, this expression can be analytically computed and 
the sensitivity correction is the same as before. A similar em- 
pirically motivat ed correction o n the e stimator expression is 
used in Eq.5 of ISchmidt et all (|2012h . where the factor is 
computed with simulations. 



2.3 Simulations 

In order to test the estimation of sizes with lensfit we 
have generated the sa me type of simulations used in 
the GREA T10 challenge ((Bridle et al.ll2010l ; IKitching et all 
I2010l . l2012ft . but with non-zero k. Multiple images were gen- 
erated, each containing 10,000 simulated galaxies in a grid of 
100x100 postage stamps of 48x48 pixels; each postage stamp 
contains one galaxy. 

Each galaxy is composed of a bulge and a disk, each 
modelled with Sersic light profiles: 



J(r) ~/oexp{ -K 



- I -1 

Td 



(6) 



where Io is the intensity at the effective radius ra that en- 
closes half of the total light and K = 2n — 0.331. The disks 
were modelled as galaxies with an exponential light pro- 
file (n — 1), and the bulges with a de Vaucouleurs pro- 
file (n = 4). Ellipticities for bulge and disk were drawn 
from a Gaussian distribution centred on zero with disper- 
sion a = 0.3. Both components were centred in the postage 
stamp with a Gaussian distribution of a — 0.5 pixels. The 
galaxy image was then created adding both components. 
The S/N was fixed for all gal axies of the image and the im- 
plementation is the same as in Kitching et al.ll|2012r h Finally 
the PSF was modeled with a Moffat profile with /3 = 3, with 
FWHM fixed for all galaxies on the image, with different 
ellipticies, drawn from a uniform distribution, with ranges 
given in Table. [Tj 

We have generated 4 different sets of images, whose de- 
tailed information is given in Table. [JJ The different types of 
image were generated to study the effects of the bulge frac- 
tion (fraction of the total flux concentrated in the bulge), the 
S/N and the PSF separately. In summary, the main charac- 
teristics of the considered sets are: 

• Set 1. Disk-only galaxies, (bulge fraction = 0), negligible 
PSF effect (FWHM PSF = 0.01 pixels) and different S/N. 

• Set 2. Disk-only galaxies, (bulge fraction = 0), with 
S/N=20 and different sizes of PSF. 

• Set 3. Negligible PSF effect, S/N=20 and different bulge 
fractions. 

• Set 4. Bulge fraction of 0.5, FWHM of PSF 1.5 times 
smaller than the characteristic size of the disk, and different 
S/N. 

To characterize the size of the galaxy the half-light disk ra- 
dius, s = Td is used. We have drawn a Gaussian distribution 
for r d , with expected value of 7 pixels and dispersion of 1.2 
pixels, to keep disk sizes of at least 2 pixels and not larger 



that the postage-stamp. The galaxy sizes explored here have 
a somewha t smal ler range (a(\nR) ~ 0.18) than found by 
IShen et all l|2003h with the SSDS catalogue, where in terms 
of pixels the mean value of the full sample is around 5 with 
aQnR) ~ 0.3 (see Fig.l of lShen et all (120031 )1. Therefore the 
sensitivity corrections are consequently larger than would be 
needed for real data. Besides a wider distribution of galax- 
ies, the important change from the original images for the 
GREAT10 challenge is the addition of a non-zero kappa- 
field that creates a size-magnification effect (in GREAT10 
only a shear field was used to distort the intrinsic galaxy im- 
ages). A Gaussian convergence field with a simple power-law 
power spectrum of Pk oc 10 _J / -11 has been applied to each 
image using the flat-sky approximation. The power-law is 
good approximation to the theor etical power s pectrum over 
the scales 10 < I < 1000 (see i.e. lHuteredl2002l) . 



3 RESULTS 

Before trying to estimate the convergence field of our most 
realistic image, we have tested the dependence on differ- 
ent aspects separately: S/N, PSF size and bulge-fraction. 
We expect these to be the observable effects that have the 
largest impact on the ability to measure the size of galaxies. 
Lower S/N will cause size estimates to become more noisy 
and possibly biased (in a similar way as for ellipticity, see 
iMelchior fc Violall2012h : a larger PSF size will act to remove 
information on galaxy size from the image and a change in 
galaxy type or bulge-fraction may cause biases because now 
two characteristic sizes are present in the images (bulge and 
disk lengths). In order to study carefully the sensitivity of 
our estimator to systematic noise, PSF or galaxy proper- 
ties, we started from the simplest case and added increasing 
levels of complexity. The number of galaxies used for the 
analysis is 200,000 for the first three sets and we increased 
the number to 500,000 for the last test to give smaller error 
bars. 

We compare the estimated k computed as in Eq.[TJwith 
the input field k, and fit a straight line to the relationship 
to estimate a multiplicative bias m and an additive bias c: 



(m + l)/t + c. 



(7) 



To obtain m and c values, we estimate the observed size of 
each galaxy, and compute k using the estimator in Eq.Q] The 
sensitivity correction is applied as explained in Sec. 12.2.11 
Then we compare the corrected estimated k of each galaxy 
with the original convergence field introduced. Finally we 
bin the datqj and a linear fit is done to compute m and c. 
This process is shown in Fig. [TJ 

We now discuss each of the categories in turn. 



3.1 Signal-to- noise 

As a first approach to the problem, disk-only galaxies with 
a negligible PSF (FHWM = 0.01 pixels) and zero ellipticity 
were generated to test the sensitivity on the S/N only given 



2 Note that the differences between the results before and after 
the binning of k are within the errorbars. 
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Set Name 


S/N 


fwhm PSF 


e PSF 


B/D fraction 




rd 








e 




Set 1 


[10,40] 


0.01 








(r d ) 


= 7 , cr = 


1.2 











Set 2 


20 


[0.01,10] 








(r d ) 


— 7,0" — 


1.2 




<e) 


= 0,cr 


= 0.3 


Set 3 


20 


0.01 





[0,0.95] 


(r d ) 


= 7 , a = 


1.2 


r-d/2 


<e> 


= 0,cr 


= 0.3 


Set 4 


[10,40] 


4.5 


[-0.1,0.1] 


0.5 


(r d ) 


= 7 , cr = 


1.2 


(r b ) = 3.5 , cr = 0.6 


<e> 


= 0,cr 


= 0.3 



Table 1. Major characteristics of the different sets used in this analysis. In bold are marked the variables explored in each set and the 
range of variation. In Set 4, PSF ellipticities are drawn from a uniform distribution in the range specified. Note that ri corresponds to 
the half-light radius. Last column is the galaxy ellipticity and is the same for both components, bulge and disk. 




'2 4 6 8 10 12 -0.1 -0.05 0.05 0.1 -0.02 -0.01 0.01 0.02 



i k. k 

n in in 

Figure 1. Sequence of steps to obtain m and c values. First panel shows the lensfit output size compared to the input size, in the second 
panel the estimated re compared to the input convergence at each galaxy, and in the third panel is shown the same plot using bins. Slope 
and intercept values of the fitting are shown in each plot (Throughout, we fit generically y = bx + c, with b = m + 1 and c = c of eq. 0. 
This is for galaxies of Set 1 with signal- to- noise 40. 



otherwise perfect data. In Fig. [5] we can see that as we in- 
crease the S/N the accuracy of the size estimation grows, as 
expected. In Fig. [3] we show the estimation of the conver- 
gence field, corrected by the sensitivity correction. Although 
the range of the output and input values differs by almost 
one order of magnitude, there is a clear correlation between 
the inputs and the outputs, and the slope is close to unity 
for all S/N explored. In Fig. [4]the estimates for m and c are 
shown, with the sensitivity correction and without it. In this 
case the correction does not alter much the results except 
for low S/N, because the sizes are less accurately estimated. 
In this paper the factor a 2 /(a 2 + b 2 ) is estimated by the in- 
verse of the slope of the size estimation fitting (see Fig. [2]). 
Using 200,000 galaxies for this test, the values found for m 
and c are consistent with zero, typically m ~ 0.02 ± 0.05, 
and c ~ (5 ± 5) x 10" 4 . 

3.2 PSF effect 

To study the uncertainties on the size estimation due to 
the PSF size, we generated images with different FWHM 
PSF values, with an intermediate signal to noise (S/N=20), 
maintaining the same properties as before, except that we 
considered here a Gaussian distribution of ellipticities with 
mean value of e = and fj e = 0.3 (per component). The size 
estimates are good for small PSFs, but become progressively 
more biased as the PSF size increases beyond the disk scale 
length (see Fig. [S]). A PSF with a FWHM larger or similar 
to the size of the disk, tends to make the galaxy look larger, 
and the estimator for k is biased. This effect can be seen 
in the slope and intercept of k vs k plot (Fig. [B]). Fig. [7] 
shows the variation of the parameters m and c with the 
ratio between the scale-length of the PSF and the galaxy 

(ratio=r d /PSFFWHivi)- 

We find no evidence for an additive bias, but we do find 
a multiplicative bias for large PSFs. With a wide size dis- 







Figure 2. Comparison of the estimated sizes by lensfit with the 
input galaxy size for different S/N in the range [10,40]. Disk- 
only circular galaxies with a negligible PSF effect are considered 
(Set 1). Slope and intercept of the fitting are shown (b and c, 
respectively). Note that the input size is the lensed one. 



tribution, some of the smaller galaxies are convolved with a 
PSF larger than its size, and this could produce an overall 
bias in k, but if the number is not very large, this will not 
affect the global estimation of the convergence field. Similar 
biases exist with shear for large PSFs, but the biases are 
larger here. For a space-based experiment, with a relatively 
bright cut at i ~ 24, such as planned for Euclid, the limita- 
tion on PSF size wil l not be a problem, as the median galax y 
size is 0.24 arcsec l|Simard et al.ll2002l : [Miller et al.ll2012T l. 
comfortably larger than the PSF FWHM of 0.1 arcsec. For 
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S/N=10 



S/N=20 



fwhm=0.0 



fwhm=1 .0 






-0.02 -0.01 




-0.02 -0.01 



Figure 3. Comparison of the binned estimated convergence and 
the input value for Set 1 with different S/N in the range [10,40]. 
Slope and intercept of the fitting are shown (b and c, respec- 
tively). For errors, see text. 







Figure 5. Sizes estimates vs input sizes for four different PSF 
scale-lengths between 0.1 and 7 pixels. Galaxies are disks with 
S/N=20 and mean size 7 pixels. Slope and intercept of the fitting 
are shown (b and c, respectively). 




10 20 30 40 

S/N 



□ c 

V c corrected 



10 20 30 40 50 

S/N 



Figure 4. m and c values computed with 200,000 galaxies of 
Set 1. Triangles are for the values obtained with the sensitivity 
correction and squares without it. 



ground-based surveys, such as CFHTLenS and future exper- 
iments the situation is not so clear, the measurement will be 
more challenging, and large empirical bias corrections of the 
order of m ~ —0.5 will be needed (see the first point of 
Fig. EI). 



3.3 Bulge fraction 

In this test we generated galaxy images with bulges with 
different fractions of the total flux, to test the response to 
the galaxy type. In Fig. [5]we show that galaxy size estimates 
for a bulge fraction of 0.2 is much better than for galaxies 
with bulge-fraction of 0.95. This is because for bulge domi- 
nated models the central part of the galaxy becomes under 
sampled due to a limiting pixel scale. The poor estimation 
of sizes is reflected in the k estimation (Fig. [5}. 

The parameters m and c for this set are shown in 
Fig. 1101 where we can see that for bulge-fraction greater 
than 0.8, the results are clearly biased (m =-0.25). For all 
bulge-fractions the error bars are around 10%. 

Although for bulge-only galaxies, the k estimates are 
poor, most of the galaxies used for weak lensing exper- 
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Figure 6. re estimates vs input values for four different PSF scale- 
lengths. Galaxies are disks with S/N=20 and mean size 7 pixels. 
Dashed line is K ou t = K; n and the solid line is the least squares fit, 
with slope and intercept shown in the plots. Note that b=m+l. 



imen ts have bulge fractions lower than 0.5 ijSchade et al.l 
1996), then the useful population can be large enough to do 
a successful analysis. 



3.4 Most Realistic Set 

The last set includes realistic values for all the effects we 
investigate. We have generated elliptical galaxies with a 
bulge fraction of 0.5 convolved with an anisotropic PSF with 
FWHM of 4.5 pixels (1.5 times smaller than the character- 
istic scale-length of the disk), and again we investigate the 
dependence on S/N. This is also a challenging test for lens&t, 
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Figure 7. m and c values computed with 200,000 galaxies of 
Set 2. Triangles represent the values obtained with the sensitivity 
correction and squares without it. 



Figure 10. m and c values computed with 200,000 galaxies of 
Set 3. Values obtained with the sensitivity correction are marked 
by triangles and by squares are without the correction. 
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Figure 8. Sizes estimates vs input sizes for different bulge frac- 
tions in the range [0.2,0.95]. Bulge+Disk elliptical galaxies are 
used, with S/N=20 and negligible effect of the PSF (Set3). Slope 
and intercept of the fitting are shown (b and c, respectively). 
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Figure 11. m and c parameters for 500,000 galaxies of Set 4. 
Squares are raw m and c values; triangles have the sensitivity 
correction included. 



the current version of which uses a simplified parameter set 
where the bulge scale length is assumed to be half the disk 
scale length. Here, we include a dispersion in the bulge scale 
length of 0.6 pixels around a mean value of 3.5 pixels. The 
analysis was done with 500,000 galaxies, to keep the error 
bars smaller than 10%. In Fig. [11] are shown the values of 
m and c for this set, with and without the sensitivity cor- 
rection. If we compare it with the previous plots we can 
see that as the galaxies get more realistic the importance 
of the correction increases. Results for this set are shown in 
Fig. 1111 showing unbiased results except for S/N=10, which 
has m = —0.19 ± 0.1. For the higher S/N points, we find 
\m\ < 0.06 with errorbars of ±0.09. 



0.03 
0.02 
0.01 

J 
-0.01 

-0.02 
-0.03 



BF=0.20 



BF=0.95 





-0.02 -0.01 





k 



-0.02 ■ — 

-0.02 -0.01 



Figure 9. k estimates vs input values for two different bulge 
fractions [0.2,0.95]. Bulge+Disk elliptical galaxies with S/N of 20 
and negligible effect of the PSF (Set 3). Dashed line is Kout = K,% n 
and the solid line is the fit of the output values. Note that b=m+l 
of cq. 



4 CONCLUSIONS AND DISCUSSIONS 

In this paper we present the first systematic investigation of 
the performance of a weak lensing shape measurement meth- 
ods' ability to estimate the magnification effect through an 
estimate of observed galaxy sizes. We performed this test by 
creating a suite of simulations, with known input values, and 
by using the most advanced shape measurement available at 
the current time, lensfd. 

A full study of the magnification effect using sizes was 
performed testing the dependence on S/N, PSF size and type 
of galaxy. The requirements on biases on shear (or equiva- 
lently convergence) for Euclid, so that systematics do not 
dominate the very small statistical errors in cosmological 
parameters are stringent: the multiplicative and additive bi- 
ases are required to be |m| ^ 2 x 10~ 3 and |c| ^ 10 -4 . A 
much larger study will be required to determine whether 
these requirements can be met for size, but we find in this 
study no evidence for additive biases at all, and no evidence 
for multiplicative bias provided that the PSF is small enough 
(<galaxy scalelength/1.5), the S/N high enough (^ 15), 
and the bulge not too dominant (bulge/disk ratio < 4). 
We find that if the disk is faint in comparison to the bulge 
(bulge/disk ratio ^ 1.5), the estimation of the half-light disk 
radius is biased. This effect increases with the bulge/disk ra- 
tio, and we find that the bias in the convergence estimation 
can be important for bulge fraction > 0.8. 

Besides instrumental and environmental issues, there 
can be astrophysical contaminants associated with weak 
lensing. In the case of shape distortion, the first assump- 
tion that galaxy pairs have no ellipticity correlation is 
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not entirely accurate. There are intrinsic alignments of 
nearby galaxies due to the alignment of angular momen- 
tum produced by tidal shear correl at ions (II correlations, se e 
for detections Brown et all (|2002|); iHevmans et alj (120041 ); 
IMandelbaum et alj (l201ll) 1; I Joachimi et al.1 d201lll2012h. an d 



for theory Heavens et al.l Jiobol); Catelan fc Porcianil (l200ih 



ICrittenden et al.l (|200lF : iHevmans fc Heavens! (120031 )). In 



addition there can be correlation s between density field s 
and ellipticities (GI cor relations, iHirata fc Seliakl (|2004h ; 
IMandelbaum et alj (120061 )). These intrinsic correlations have 
been studied in detail and it is not trivial to account for or 
remove them when quantifying the weak lensing signal. 

The intrinsic correlation of sizes and its dependence on 
the environment, are still open issues. In fact, the corre- 
lation of sizes and density field, is known to play an im- 
portant role in discriminating between models of size evolu- 
tion; recent work finds a significant correlation between sizes 
and the density field using around 11, 000 galaxies drawn 
from the joint DEEP2 /DEEP3 data-set (|Cooper et al.ll2012l ; 
IPapovich et~aH |2012| ), while earlier studies with smaller 
samples have been in disag reement. Using 5,000 galax- 
ies of STAGES data-set, iMaltbv et ail (|2010l ) find a pos- 
sible anti-correlation between density field and size for 
intermediate/low-mass spiral galaxies. Clustered galaxies 
seem to be 15% smaller than the ones in the field, while they 
do not find any correlation for high-mass galaxies. Also for 
mas sive elliptical galaxie s from ESO Distant Clusters Sur- 
vey, iRettura et al.l ([2010) do not find any significant corre- 
lation, while using the same data set Cimatti et alj (120121 ) 
claims a sim i lar co rrelation as in ICooper et al. I I2OI2I ). In 
iPark fc Choil (|2009h they study the correlation between sizes 
and separation with late and early-type galaxies from SSDS 
catalogue, at small and large scales. They compare the size 
of the nearest neighbour with the separation between them, 
and find larger galaxies at smaller separations. This correla- 
tion is found for early-type galaxies if the separation between 
the galaxies is smaller than the merging scale, but not for 
larger separations. The size of late-type galaxies seems not 
to have a correlation with the separation in any scale. We 
expect that further studies with larger samples will clarify 
the intrinsic correlations of sizes, we note that the system- 
atics are generated from different physical processes than in 
the case of shear and so will affect the signal in a different 
way; we suggest this is a positive, and another reason why 
a joint analysis of ellipticity and sizes is interesting. 

The analysis presented here has assumed that the statis- 
tical distribution of galaxy sizes is known, whereas in prac- 
tice the size distribution depends on galaxy brightness and 
must be determined from observation. Gravitational lensing 
of a galaxy with amplification A increases both the inte- 
grated flux and area of that galaxy by A, which has the 
effect of moving galaxies along a locus of slope 0.5 in the 
relation between log(size) and log(flux). Thus, if the intrin- 
sic distribution of sizes r of galaxies scales with flux S as 
r oc 5 1 ' 9 , the apparent shift in size caused by lensing amplifi- 
cation A is r' oc A°- 5 ~P ; resulting in a dilution of the signal 
compared with the idealised case investigated in this pa- 
per. A similar effect occurs in galaxy number magnification, 
where the observed enhancement in galaxy number density 
N' varies as N' oc A 01-1 , if t he intrinsic number den sity of 
galaxies varies as N oc S~ a jBroadhurst et al.| [T995). The 
value of /3 at faint magnitudes has recently been estimated 



bv lMiller et al. J2012l). who a nalysed the fits to galaxies with 
i < 25 of lSimard et alj (|2002l ) and estimated /3 ~ 0.29. Thus 
we expect this effect in a real survey to dilute the lensing 
magnification signal by a factor 0.42, but still allowing detec- 
tion of lensing magnification. In practice, the dilution factor 
could be evaluated by fitting to the size-flux relation in the 
lensing survey. 

Lensing number magnification surveys are also affected 
by the problem that varying Galactic or extragalactic ex- 
tinction reduces the flux of galaxies and thus may cause a 
spurious sig nal (e.g. iMenard et al.l l20ldh . Such extinction 
would also affect the size magnification of galaxies, but with 
a different sign in its effect. Thus a combination of lensing 
number magnification and size magnification might be very 
effective at removing the effects of extinction from magnifi- 
cation analyses. 

Space-based surveys as Euclid should overcome the lim- 
itations that we have exposed here, having a large number of 
galaxies, with S/N > 10, and importantly a PSF at least 1.5 
smaller than the average disk size. The addition of the size 
information to the ellipticity analysis is expected to reduce 
the uncertainties in the estimation of weak lensing signal, 
and therefore improve the constraints of the distribution of 
matter and dark energy properties. 
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